Nature of the spin dynamics and 1/3 magnetization plateau in azurite 
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We present a specific heat and inelastic neutron scattering study in magnetic fields up into the 
1/3 magnetization plateau phase of the diamond chain compound azurite Cu3(C03)2(OH)2. We 
establish that the magnetization plateau is a dimer-monomer state, i.e., consisting of a chain of 
5 = 1/2 monomers, which are separated by S — dimers on the diamond chain backbone. The 
effective spin couplings Jmono/ks = 10.1(2) K and Jdimer/kB = 1-8(1) K are derived from the 
monomer and dimer dispersions. They are associated to microscopic couplings Ji/fcs = 1(2) K, 
Ji/kB = 55(5) K and a ferromagnetic J^/ks = —20(5) K, possibly as result of d z 2 orbitals in the 
Cu-0 bonds providing the superexchange pathways. 

PACS numbers: 75.30.Et, 75.10.Pq, 75.45.+j 



Great interest has surrounded the observation of a 1/3 
magnetization plateau in azurite Cu 3 (C0 3 ) 2 (OH) 2 [E0- 
This material, famous as a painting pigment of deep- 
blue colour, has been proposed as a realisation of the 
exotic diamond-chain Hamiltonian of coupled spin-1/2 
moments, written as 
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Here, J 2 is the magnetic coupling of the diamond back- 
bone, while Ji and J 3 represent the coupling of the 
monomers along the chain 

BUS (Fig. E|). Depending 
on the relative coupling strengths Ji , J%, J3, this model 
affords a host of exotic phases and quantum phase transi- 
tions, including possibly M = 1/3 fractionalisation Q or 
exotic dimer phases [4]. However, determining the mag- 
netic exchange couplings in azurite has proved difficult, 
yielding controversial results. While a susceptibility x 
study claims Ji/k B = 19 K, J 2 /k B = 24 K, J 3 /k B = 8.6 
K, implying strong frustration [l[, subsequent numerical 
studies of \ dispute this claim, proposing a ferromagnetic 
(FM) J 3 , and thus a non-frustrated scenario Q. 

The general issue underlying these starkly contrast- 
ing interpretations of the same experimental data is that 
of the nature of magnetic coupling in low-dimensional 
(low-D) quantum magnets. In azurite Cu 3 (C0 3 )2(OH) 2 , 
the Cu 2+ ions [S = 1/2) are in a square-planar coor- 
dination on two inequivalent sites Q. The system has 
a monoclinic crystal structure (space group P2i / c, lat- 
tice parameters a = 5.01 A, b = 5.85 A, c = 10.3 A, 
(3 = 92 A [8, !|), where the Cu 2+ network is built up by 
diamond-shaped units arranged in chains running along 
the b direction. The magnetic exchange pathways for the 
Ji, J 2 and J 3 couplings are along Cu-O-Cu bonds with 



angles of 113.7°, 97° and 113.4°, respectively. Hence, 
the Goodenough-Kanamori- Anderson rules for superex- 
change [lOj would predict the couplings to be weakly an- 
tiferromagnetic (AFM) [111 ]. 

However, this conclusion is only valid if the magnetic 
orbitals are dominantly of d x -i_ y i character. If the Cu or- 
bitals tend to be more d z i like, higher-order effects come 
into play, which can give rise to FM exchange [l2]. The 
latter case may become relevant for systems with bond 
angles away from the limiting cases of 90° and 180°, 
such as CuO or some molecular magnets [l3j. In fact, 
recent x-ray diffraction experiments [9j on azurite indi- 
cate a dominant d z i character of Cu 2+ , making azurite 
another prominent candidate for such superexchange in- 
teractions. 

In this Letter we present a detailed study of azurite by 
means of specific heat and inelastic neutron scattering 
in zero and applied magnetic field. We investigate the 
elementary magnetic excitations, in particular the q de- 
pendence of the spin-excitation spectrum in the plateau 
phase, and estimate the microscopic coupling constants. 

Specific heat C P {T) measurements at temperatures 1.6 
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to 30 K were carried out using an ac-calorimeter 
an azurite crystal (mass: 0.36 mg), which was cut from 
the crystal used for neutron scattering. An external field 
up to 4 T was oriented in the ac-plane, and 65° away from 
the c axis. This orientation corresponds approximately 
to the easy axis of the AFM phase below Tn = 1.85 K 
0. The C P (T) results are displayed in Fig. Q] up to 10 
K. Beside the AFM transition, indicated by a sharp dis- 
continuity, the zero-field data reveal a broad maximum 
around 3.7 K. At B = 4 T the maximum becomes re- 
duced in size and shifted to lower temperatures. 

The starting point of our analysis is the dimer- 
monomer model proposed in Ref . [4] , with the intradimcr 
coupling constant J 2 representing the dominant energy 
scale. At temperatures T < 10 K, considered here for 
the specific heat, and for J 2 //cs 3> 10 K (see below), the 
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FIG. 1: (Color online) C P (T) of a single crystal of azurite 
for B = T (squares) and 4 T (Q)- The solid (dashed) 
line represents a fit to the zero-field (B — 4 T) data using a 
magnetic (C ma3 ) and phononic (Ci at ) contribution. The 
inset shows C m ag{T) at B = T, C; a t and the magnetic 
entropy S mag derived by integrating (C P (T) ~ Ci a t(T))/T. 

magnetic degrees of freedom of azurite are dominated by 
the chain of spin-1/2 Cu 2+ -monomers, which are anti- 
ferromagnetically coupled via the rungs of the diamond 
backbone [lj. Using the magnetic specific heat C mag of 
the AFM S =1/2 Heisenberg chain (AFHC) and in- 
cluding a lattice contribution Ci a t <x (T/0c) 3 , the zero- 
field specific heat C P (T) is well fitted down to 2.5 K. In 
this fit, the magnetic coupling Jafhc I^b = 7.0(1) K and 
a Debye temperature Qd = 188 K were used (solid line in 
Fig. [T|). The corresponding magnetic entropy S mag , nor- 
malized to a Cu 2+ monomer, reaches the full value Rln2 
at about 10 K (broken line in the inset of Fig. [T]) . Us- 
ing the same lattice contribution and C maa {T, B) results 
calculated for the AFHC in finite field [161], this ansatz 
(broken line in Fig. QJ also captures the main features of 
the experimental data at B = 4 T. 

Low-temperature (T = 1.5 K) inelastic neutron scat- 
tering (INS) was used to further probe the character of 
the magnetic excitations both at zero and high magnetic 
field. Experiments were carried out at the Berlin Neutron 
Scattering Center, BENSC, using the V2 cold-neutron 
triple-axis spectrometer. Constant-Q, energy transfer 
scans in a range - 7.8 meV were performed with inci- 
dent neutrons fixed at ki = 1.55 and 1.3 A" 1 which gave 
an energy resolution of 0.15 and 0.11 meV, respectively. 
A large (> 15 g), naturally grown azurite single crystal 
was mounted with a horizontal a* — b* scattering plane, 
and magnetic fields of up to 14 T applied perpendicular 
to b*. Susceptibility data, taken on pieces of this crystal, 
perfectly reproduce the results reported in Ref. [l|. 

The essential features of our INS studies are summa- 
rized in Fig. [2] Here, we plot the energy dependence of 
the scattering intensity for various positions q = [1 k 0] 
along the diamond chain in zero field and in 14 T, i.e., in 
the 1/3 magnetization plateau phase between B c \ = 11 
T and B C 2 = 30 T. Data shown here only cover the range 



< k < 0.5 since measurements up to k = 1 reveal 
symmetric behaviour with respect to k = 0.5. 
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FIG. 2: (Color online) Energy dependence of the INS spectra 
of azurite at q = [1 k 0], k = — 0.5, in zero field (left) and a 
field of 14 T 1 6* (right) at T = 1.5 K. Lines indicate fits to 
the data, which are shifted for clarity; for details see text. 

Along q = [1 k 0] and in zero magnetic field, there 
is scattering intensity at energies < 2 meV displaying 
the | sinq| dependence of an S — 1/2 AFHC (left panel 
Fig. [2]). The fits shown in the figure were produced by 
combining Gaussian line shapes for the features above 2 
meV, while those below 2 meV were calculated using a 1- 
D chain [13] convolved with instrumental resolution and 
the copper form factor. The fit yields an effective mag- 
netic coupling strength of Jafhc I^b — 9(1) K (solid 
lines < 2 meV). The difference between Jafhc measured 
from specific heat and inelastic neutron scattering is not 
well understood and will require additional studies. Per- 
pendicular to the chain, i.e., along [h 0.5 0], within ex- 
perimental resolution there is no dispersion (not shown) , 
consistent with the 1-D nature of these spin excitations. 

Further, our spectra reveal a broad peak structure at 
energies above 3.5 meV in zero magnetic field (left panel 
Fig. [2|) . As we will argue below, this feature arises from 
singlet-triplet excitations of the spin dimers on the back- 
bone of the diamond chain structure. 

Within the dimer-monomer model, the monomer chain 
becomes fully polarized at the lower critical field B c \ of 
the plateau phase. Correspondingly, a gap opens at the 
AFM point q = [1 0.5 0]. The excitations are now fer- 
romagnons, yielding scattering intensity with a cosinuso- 
dial dispersion (see low-energy peak 0.4 < E < 2.2 meV 
in 14 T; right panel Fig. [2]). By fitting the peak posi- 
tions using Gaussian profiles as indicated in the plot, we 
extract the E{k) dependence in Fig. [3] (black squares). 
The data are parameterized by an expression E(k) = 
g av ^BB + J mono {cos(2nk) - 1) + A cham . With g av = 2.1 
from ESR measurements (l8j and A c h a in = 0.53 meV, 
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we obtain the effective magnetic coupling Jmono/ks — 
lO.f (2) K. Using higher order cosine-terms does not im- 
prove the fitting, resulting in 0.2 K as upper boundary 
for next-nearest-neighbor interactions. 
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FIG. 3: (Color online) (top) Dimer-monomer model with mi- 
croscopic couplings Ji, J2, J3, and effective couplings Jmono, 
Jdimer, characterizing the plateau phase, (bottom) E(k) de- 
pendence along q — [1 k 0] of azurite for the two low-lying 
excitations in B(_L &*) = 14T. Lines are fits to the data, see 
text. The small asymmetry of the data with respect to k = 0.5 
reflects a slight misalignment of the crystal. Highlights I, II, 
III indicate positions investigated for their B dependence in 
Fig.H 

The broad peak in the spectra above 3.5 meV in Fig. 
[2] is attributed to singlet-triplet excitations within the 
dimers. In magnetic fields, Zeeman splitting lifts the 
zero-field degeneracy of the singlet-triplet peak, lowering 
the Iff) triplet branch in energy and separating it from 
the broad structure. Accordingly, in B = 14 T, a weakly 
dispersive peak is observed at ~ 2.2 — 2.5 meV (right 
panel Fig. [2|) . By fitting the peak positions of the low- 
lying triplet using Gaussian profiles, we extract the E{k) 
dependence in Fig. [3] (red circles). It can be parameter- 
ized by E(k) = -QavUBB + J 2 + JdimerCOs(2Tlk) + /J B 6, 

with g av =2.1 [18[ and J2 as zero-field singlet-triplet 
splitting. The parameter b accounts for the internal 
field shift at the dimer site due to the alignment of 
nearest-neighbouring monomers. From the width of the 
dispersion we obtain an effective dimer-dimer coupling 
Jdimer = 1-8 K. The other triplet branches, even in high- 
est fields, superimpose and yield a broad distribution of 
scattering intensity for energies > 3.5 meV, thus pro- 
hibiting a precise determination of peak positions. As 
yet, the cause for the peak broadening of these branches 
is not understood. 

The critical fields of the plateau phase are extracted 
from a field-dependent study of the monomer and lowest 
dimer excitations. First, we note that the field depen- 
dence of the monomer dispersion should be that of the 
5 = 1/2 AFHC. Hence, at q = [1 0] (marker I in Fig.EJ) 
a (close to) linear-in-field behaviour is expected, which is 
gapless in zero field [19j. As shown in the lower left panel 



of Fig. [4] a linear behaviour is in fact observed experi- 
mentally. Further, B c \ = 11 T is identified in the field 
dependence of the peak intensity (upper left panel Fig. 2]) 
as a distinct kink. The intensity is constant above B c i, 
but it falls off rapidly with decreasing field for B < B c \. 
The intensity decrease is due to the first moment sum 
rule and is related to the magnetization [20| . 
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FIG. 4: (Color online) (left) B dependence of the integrated 
peak intensity (top) and monomer energy (bottom) at q = 
[1 0] (I) and [1 0.5 0] (II). (right) Energy of the low-lying 
triplet branch C I T T> ) at q = [1 0.5 0] (III) as function of B. 
Hatched region indicates the plateau phase. Solid lines: fits 
to the data; dashed lines: guides to the eye. 

Surprisingly, an extrapolation of the peak position for 
B > B c i to zero field yields a finite gap A c h a in = 0.53 
meV. This feature also accounts for the mismatch at low 
energies between fit and data at q = [1 0] in zero field 
(left panel Fig. [2]). Unfortunately, a detailed study of this 
behaviour is hampered by the strong reduction of peak 
intensity (upper left panel Fig.rjJ. Possibly, the AFM or- 
dered state affects the behaviour of the monomer chain, 
leading to a residual zero-field gap at the nuclear zone 
centre. Certainly, this issue calls for further investiga- 
tions. 

In addition, we measured the B dependence of the 
monomer peak at [1 0.4 0]. Assuming a cos(2irk) de- 
pendence of E(k), we approximate the field dependence 
of the E(k) minimum at [1 0.5 0] (marker II in Fig. [3]), 
this way accessing the field of gap closure for the ferro- 
magnons [2l[. For B > B c \ the expected linear-in- field 
behaviour of the gap is observed experimentally, with an 
extrapolated field of gap closure Bo = 9.7 T (lower left 
panel Fig. [4|. We associate the difference between B c \ 
and Bq to the local mean field due to the AFM order. 

Our data suggest that the magnetization plateau 
reaches from the field B c \ of full monomer polarization 
up to B C 2, where the low-lying triplet |ft) reaches zero 
energy. Moreover, Zeeman splitting results in a depen- 
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dence oc B for this triplet branch. To test this, we have 
measured the field dependence of the ff) branch at the 
minimum q = [1 0.5 0] (marker III in Fig. [3J data in 
Fig. [3J right panel). As expected, for fields B > B cl a 
linear field dependence is observed. An extrapolation to 
zero energy from the low-lying triplet Iff) at q = [10. 5 0] 
yields a value of 31 T, which closely matches the reported 
value B C 2 = 30 T as upper critical field of the magnetiza- 
tion plateau. Extrapolating to zero energy for q = [1 0] 
gives B C 3 = 33.6 T also comparable to the reported value 
m of 32.5 T for the onset of magnetization saturation to 
the full moment. 

Below B c i, for the triplet |ff) at q = [1 0.5 0] again a 
behaviour oc B, but with a steeper negative slope, is ob- 
served. This reflects the presence of local molecular fields 
as a result of the AFM aligned monomers. Fitting the 
data at B < B ci yields a zero field singlet-triplet split- 
ting of 4.8(5) meV (~ 55(5) K). This splitting essentially 
corresponds to the magnetic coupling Ji on the diamond 
backbone. 

Altogether, our data reveal that the magnetization 
plateau in azurite arises out of the dimer-monomer state 
proposed in Ref. Q and indicated in the top panel of 
Fig. [3] There is a chain of antiferromagnetically cou- 
pled S = 1/2 monomers, which are fully polarized in the 
plateau phase. The coupling between the monomers is 
provided by the rungs in the diamond units, which carry 
two spins forming a dimer (Ji)- Nearest- neighbor dimcrs 
are weakly coupled (Jdimer) to each other. 

Within first-order perturbation theory, analytical ex- 
pressions relate Ji, Jdimer to the microscopic couplings 
Ji , J3 , provided that J2 S> J\ , J3 [5] . Although in azu- 
rite the latter condition might not be strictly fulfilled, 
our data indicate J2 to be the strongest coupling. Em- 
ploying J2 and Jdimer from experiment, we approximate 
Ji and J3 following Ref. [H. First, for the dimer cou- 
pling it predicts Jdimer ~ {Ji — J3) 2 I^Ji- Secondly, the 
change in slope of the dimer dispersion at B c \ (Fig. [4|) 
arises from a change in local mean fields from AFM to 
FM. Extrapolating the data at B > B c \ to zero field 
we obtain the difference between FM and AFM polariza- 
tion AE k (Ji + J 3 )/2 + Jdimer = -0.8 mcV. Finally, 
we can use the various critical field expressions such as 
B C 2 ~ {J2- Jdimer + {Ji + J3)/2) /gav^B for checking the 
consistency of our calculations. With this procedure we 
obtain J 3 /k B ~ -20(5) K and J x /k B ~ 1(2) K. 

Most importantly, our findings indicate that there is 
a FM coupling in azurite, as was suggested previously 
Q. The J3 path has Cu-O-Cu exchange along one filled 
and one unfilled orthogonal orbital, with atomic Cu-0 
distances of 1.947 and 1.989 A separated by an angle of 
113.4°. By reducing the Anderson superexchange mecha- 
nism and with the Cu orbitals of a dominantly d z 2 charac- 
ter Q , direct exchange and oxygen Hund's rule coupling 
can yield a substantial ferromagnetic exchange pathway 
[I3 |. Further, with |Ji| <C \J2\, \Js\ frustration effects will 



be absent or weak. J\ represents a weak perturbation, 
and is possibly less important than other factors such as 
Dzyaloshinskii-Moriya-interaction []J . 

To conclude, we have studied the distorted diamond 
chain model system azurite CU3 (003)2 (OH)2 by means 
of specific heat and inelastic neutron scattering. The 
magnetization plateau can be understood within the 
dimer-monomer model from Ref. |4J . For the microscopic 
couplings Ji, Ji and J3 within the diamond units we 
find Ji and J3 to be dominant, with J3 being ferromag- 
netic. Accounting for this type of superexchange ought 
to be possible along the lines set out in Ref. [l2|. Given 
the limits of our perturbative analytic approach, fur- 
ther calculations, which take the renormalisation of mag- 
netic couplings by quantum fluctuations into account, are 
needed to fully model the distorted diamond chain azu- 
rite Cu 3 (C0 3 )2(OH)2 0- 
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